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Davey-Stewar tson  I  is  a  nonlinear  evolution  equation 
originally  derived  in  the  context  of  multidimensional  weakly 
nonlinear  water  waves.  It  has  recently  been  exactly  solved  by 
the  classical  inverse  scattering  method  for  localized 
potentials,  and  also  possesses  nonlocal  soliton  solutions.  We 
have  calculated  Poisson  bracket  relations  for  elements  of  the 
scattering  matrix  ,  as  well  as  corresponding  quantum  commutation 
relations.  Commutation  relations  are  found  that  are  a  2+ld 
generalization  of  a  Yang-Baxter  algebra.  > 
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Exactly  solvable  systems  have  played  a  significant  cole  in  our 

understanding  of  nonperturbat ive  phenomena  in  physics.  Many 

quantum  field  theories  in  1+1-dimensions  have  been  found  to  be 

integrable,  enabling  the  calculation  of  exact  S-matrices  and 

physical  spectra.  The  Ising  model  and  other  exactly  solvable 

models  of  2-dimensional  statistical  mechanics  have  helped  to 

provide  a  basis  for  modern  scaling  theory.  Moreover,  some  of  the 

more  interesting  mathematics  occurring  in  quantum  string  theories, 

incuding  loop  spaces  and  Kac-Moody-Virasor o  algebras,  also  appear 

% 

in  integrable  systems. 

Associated  with  every  known  integrable  quantum  system  in  2 

dimensions  (or  equivalently  1  space  +  1  time  dimension)  is  a 

1  2 

solution  of  the  so-called  Yang-Baxter ( YB )  equations  '  ,  and  the 

existence  of  a  corresponding  "Yang-Baxter  algebra”.  The  YB 
equations  arise  in  various  contexts,  and  have  come  to  be  regarded 
as  the  criterion  for  exact  quantum  integrabi 1 i ty . 

It  is  certainly  of  interest  to  extend  the  study  of  quantum 
integrability  to  higher  dimensions.  Progress  was  made  in  this 
direction  when  Zamolodchikov  considered  the  scattering  of 
"Straight  strings”  in  a  plane,  and  wrote  down  a  3d  generalization 
of  the  YB  equations,  called  the  tetrahedran  equations,  as  well  as _ 
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a  con 


jectured  solution"5.  Baxter  was  able  to  verify  this  solution  ^ 


and  exactly  calculate  the  free  energy  of  an  equivalent  classical  q 


statistical  mechanical  model  However,  the  physical 
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interpretation  of  the  Zamolodchikov-Baxter  solution  is  somewhat 
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problematic,  and  the  tetrahedran  equations  are  so  complex  that 
little  progress  has  been  made  in  finding  other  solutions. 

An  alternative  approach  to  searching  for  new  quantum 
integrable  systems  in  higher  dimensions  is  to  exploit  our 
knowledge  of  existing  integrable  classical  systems.  There  now 
exist  a  number  of  nonlinear  evolution  equations  in  (2+l)d  which 
are  solvable5  by  the  classical  inverse  scattering  transform 
(CIST)^.  According  to  our  experience  in  1+ld,  each  of  these 
classical  systems  should  have  a  corresponding  quantum  analog 
which  is  exactly  integrable.  In  this  letter  we  consider  the 
quantum  analog  of  such  a  classical  system,  known  as 
Davey-Stewartson . 
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The  Davey-Stewar tson { DS )  equation  is  a  nonlinear  partial 

differential  equation  in  2  +  ld,  originally  formulated  as  a  model 

to  describe  the  evolution  of  mul t i -d imens i ona 1  weakly  nonlinear 

water  waves.  Depending  on  the  choice  of  the  parameters  in  the 

equation,  it  admits  two  types  of  soliton  solutions,  localized 

lump-like  solitons  and  nonlocalized  straight  line-like  solitons 

Classically,  the  asymptoptic  scattering  of  the  lump  solitons  is 

trivial,  but  the  line  solitons  experience  a  nontrivial  phase 

shift.  For  our  purposes  DS  is  an  obvious  choice  because  it  is 

one  of  the  simplist  of  the  known  higher  dimensional  integrable 

classical  systems.  It  reduces  in  the  1+ld  limit  to  the 

well-known  nonlinear  Schr od i nger ( NLS )  system,  whose  quantum 
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version,  the  de 1 ta - f unct i on  gas  model  ,  or  quantum  NLS  model 


is  one  of  the  best  understood  of  the  quantum  integrable  systems. 

We  have  calculated  various  classical  Poisson  bracket  relations 
between  elements  of  the  scattering  matrix  of  the  underlying 
linear  problem  for  DSI#  enabling  us  to  explicitly  identify  the 
classical  action-angle  variables.  We  have  also  calculated 
certain  corresponding  quantum  commutation  relations  and  find  them 
to  be  a  2+ld  generalization  of  a  Yang-Baxter  algebra. 


We  first  discuss  the  classical  case.  We  will  be  concerned  with  the 
hyperbolic  version  of  the  DS  equation,  a  non-linear  partial  differential 
equation  for  a  complex-valued  function  q  •-  q(x,y,t), 

2  2 

1  ft  =  ‘  +  +  lAi<l  -  1clAo * 

L  L  3x^  3y^  1  * 

where  -^Ai  =t(37  + 


^  3x  +  3y^A2  =  2^3X  '  3y^rq^’ 


with  r  *  ±q  (q*  denoting  the  complex  conjugate  of  q). 

This  time  evolution  equation  for  q  can  be  generated  by  a  non-local 
Hamiltonian  (which  will  depend  on  the  choice  made  for  A.  and  A^)  via  the 
Hamiltonian  formulation  of  classical  mechanics,  where  q  and  r  are  the 
conjugate  variables. 

As  is  the  case  for  all  nonlinear  POE'S  solvable  by  the  CIST,  (I)  appears 
as  the  compatibility  condition  for  two  underlying  linear  equations. 


4  v  =  A  ^  +  iQ  ~ip  +  ij  —j  ip  , 


-here  Q  =  (“  J)  ■  J  ■  (J  .?) 


?<qx  *  y 


( r  -r  ) 
.  2  x  y 


and  <p  =  ip(x,y,t)  is  a  2x2  solution  matrix. 

The  first  of  these  equations,  (3a),  can  be  viewed  simply  as  a  linear 
scattering  problem  in  which  q  plays  the  role  of  the  potential.  (3a) 
for  suitable  choice  of  boundary  conditions,  can  oe  rewritten  as  a  system 
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of  linear  integral  equations, 

2i(<c„+A)J.c.  r  I 

*rU,  <,  a)  =  5  ■  j  e  R  J  J  +  j  <U'  G^U-s',  <)(QU  )*  U\  <,  a  )  )i  j  , 


where  ^  =  (x+y)  and  e2  =  (x-y)  w1th  ^  denotin9  the  coordinate  pair  ?2), 

*  =  +  iicj  is  a  complex  parameter,  A  is  a  real  parameter,  the  indices  i,  j 

can  each  take  on  values  1  or  2  (where  we  use  the  notation  1=2  and  2  ^  1,) 
and  all  integrations  are  over  infinite  space.  Also,  for  convenience  of 
notation  we  use  *»  t)  =  A,  tje^R  ^  and  we  shall 

suppress  the  argument,  t. 

We  choose  the  Greens  function 


Gijl:- <]  ~~  Gi(t'  ;ij> 


-  R  T  R 

with  <..=<..+  i<] .  <  .  =  < 

1J  ij  lj’  ij 


I  +  J1-Jj(,fR-<i  and  <■  j  -  <j 


Gu,  <)  is  obtained  by  taking  the  appropriate  limit  of  the  Greens  function 
of  the  more  general  D-bar  problem.^ 

We  also  will  find  it  useful  to  define  a  solution,  of  an  adjoint  linear 


problem. 


t  ,  L  (  £  »  L.  *- 


’i k'  *  kj 


-2i(<R  +A’)j  r  ff 

KJ  +  J  I  d?: 


^  -1  'S'  '  <g,A  )QHk^  ^Gk(t’  ‘kj 


Of  fundamental  interest  in  both  the  classical  and  the  quantum  problem  is 
the  "scattering  matrix"  or  the  "scattering  data"  of  (6),  which  we  define  to  be 
rr  "2i  (k?.  +  A'  )J  .  e. 

Ti j ( < ,  A,  x')  =  jj  dee  J  (Q(0  *U,  <,  a) ) i j  . 


For  certain  choices  of  the  parameters  <c,  a  and  A1,  T  can  be  shown  to  have 
a  very  simple  time  dependence,  and  is  thus  used  in  the  CIST  to  "reconstruct1 
the  potential  q(x,y,t)  at  arbitrary  times,  for  appropriately  given  initial 
conditions. 

We  can  calculate  Poisson  bracket  relations  between  elements  of  T,  where 
we  define  canonical  Poisson  brackets 


if,  g}=  i 


A.  6f  6q  6f 

d5  6q(  C )  6r(l)  '  4rU) 
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We  find,  by  use  of  the  linear  integral  equations,  (6)  and  (8),  that 

A,  X  ),  T^  (  t  ,  u,  u  )} 

2  rf 

_  d?  caa^’  ^ae ’  X'^  ^ae^’  <’  xa6’  *a<5^’  T’ 
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The  solution  p  and  its  adjoint  c  satisfy  Z  —  CikU.  <kj»  A1) 

k  ■*  1  k 

i'kjU.  T,  u)  =  0. 


This  identity  can  be  used  to  rewrite  the  integrand  appearing  in  (  ll) 
as  follows: 


{Ta6(<,  A'),  T^(t.  U,  u') 


-  Z 

a  =  l 


dc. 


d^  0<Ja(V5a,))  ^a(V  <ae’  x'}  *aB(^*  Ci*  X) 


Sa^a’  ^a ’  Ta6’  u'  ^  *a6^a’  Ca’  T> 


2 

+  Z  J. 


a=l 


dca  Sxa^*  <a8*  x'  ^  *a6^’  T> 


JiCS  = 


dca  Sa(5'  ’  *a6’  J '  }^a8 ( ^ '  *  <*  x) 


v;  =  - 


In  order  to  evaluate  (12)  it  is  necessary  to  find  asymptotic  expressions 


for  i*/  and  However,  these  can  be  found  easily  by  using  (6)  and  (8)  and 
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noting  that  it  is  possible  to  write  6  in  the  two  alternative  forms  (7b)  or 
(7c).  Then 

2i  (<d  +  A)J, 


lim  'K.-U,  <,  x)  =  5.  .e 

H  - 


'R  '  ksk 


2i(^Vi)J1e1 


Ho  -  «-  i\*.  •  ■  s. 

g  3(+Jki)e  kj  T,.(<,  x,  0 
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kj 


(  13) 


and 


lim  CikU,<kj,  X')  -  5. k e 
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Inserting  (13)  and  (14)  into  (12),  we  arrive  at  an  expression  for 
(T  .(<,  X,  x'),  T  .(t,  u,  u')}purely  in  terms  of  T's. 

do  yo 

Instead  of  writing  down  the  a  lengthy  expression,  which  contains  terms  up  to 
quartic  in  T,  we  instead  give  results  in  two  interesting  limiting  cases. 

First,  letting  X=x‘=u=u'=0  and  T(<,  0,  0)  =  TL(<),  we  recover  the 
scattering  data  of  the  hyperbolic  limit  of  the  D-bar  problem,  and  making 
use  of  an  identity  easily  derived  from  (6),  find  Poisson  bracket  relations 


(Tl[1(<),  T^(t)} 

(T^U),  T^(t)} 
(T^(<),  T^(t)} 


{T^2(<) ,  T^(t)}  =  {T^j(k),  T^(t)}  =0 


(  2tT  )  <5  (  J  *T  J  )  ) 


i—  +  2n  9(<i-ti) 


<r-Ti2-1£ 


T12(t) , 


as  well  as  a  number  of  other  similar  relations. 

Alternatively,  we  can  take  the  limit  <j  -*•  +  °°,  <R  -  +  ”, 

T(<,  X,  x'  )  -  T+(G,  0* ),  where  0  =  <R  +  X,  0'  =  <  R  +  x*  are  kept  finite. 

In  this  way,  we  recover  the  scattering  data  associated  with  a  solution  to 

-  2i  ©  J  .  ?  .  +  i02J-t 

(3a),  u+-U,  0)  =  4-ii(C,©)e  J  J  analytic  in  the  upper-half 

0  plane,  which  is  used  in  the  Riemann-Hi 1 bert  approach  to  CIST.  We  find 


ts*6(0,  3'),  St,  (♦.  ♦'))  =  ±s*„U,  e'l  SI5(0-  ♦')  •  (J 


y6 


2  Y  8 
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(15a) 

(15b) 

(15c) 
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4  J 


Sox  3, 


do 

s+ 

2ni  (o+i  £ ) 

a8 

(  d  o 

S  + 

2iri  (o  +  ie ) 

y6 

s+  (e  -  o,  4>' ) 

y  o 

s!q(o  ,  <t>'  +cr) , 


(16) 


Vv: 


.  A  S'. . 
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where  we  '  ve  defined  S*  (  0  ,  ©  '  )  r  2  tt  J  ${  9  -  6  )  +  T*  (©,©’). 

Similarly,  the  limit  k*  -* -«o  ,  k*-*-*-*>  gives  us  the  scattering 
data  associated  with  a  solution  analytic  in  the  lower-half 
plane . 

From  the  Poisson  bracket  relations  for  the  elements  of  TL  , 

L  l 

given  by  (15),  one  can  see  that  T^2  (  K  )  and  T21  (  K  u  )  are 

related  to  the  original  conjugate  variables,  q  and  r,  by  a 
canonical  transformation.  Furthermore,  T^L  (K  )  and  T22L  (K  ) 
are  action  variables  of  the  theory,  and  it  can  be  shown  by  using 
the  appropriate  limit  of  the  D-bar  formalism  that  they  can  both 
be  expanded  in  powers  of  l/kR  /  generating  two  infinite  sets  of 
constants  of  the  motion,  including  as  members,  the  total  momentum 
and  energy-Hami ltonian  of  the  system.  Thus  the  nonlinear 
evolution  equation,  (1),  (2),  is  exactly  integrable,  with  proper 

boundary  conditions  (D-bar). 

The  results,  (16),  for  the  Riemann-Hi lbert  formulation  of  the 
inverse  problem  have  a  very  different  form  and  reduce  in  the  1+ld 
limit  to  well-known  classical  relations  of  the  YB  type.  This 
leads  us  to  search  for  analogous  quantum  commutation  relations. 

We  have  been  able,  in  fact,  to  calculate  such  relations  for  an 
operator  version  of  DSI  with  ordering  taken  to  be  as  it  appears 
in  ( 1 ) — ( 3 ) .  (Note  that  we  do  not  treat  the  normal-ordered  case.) 
If  we  replace  the  canonical  Poisson  bracket  relations  for  q  and  r 
by  their  corresponding  canonical  equal-time  commutation 
relations,  we  can  modify  our  calculation  slightly,  now  taking 


care  to  preserve  correct  operator  ordering  throughout,  and  derive 
commutation  relations  for  elements  of  S+  (or  similarly  for  S  ). 
We  find  generalized  Yang-Baxter  commutation  relations  for  S*' 
given  by  (16),  with  (S+  ,S  +  >  replaced  by  (S+  ,S  +  )  on  the 
left-hand  side  of  the  equation. 

In  the  quantum  case  as  well  as  the  classical  case,  it  is  not 
obvious  how  to  identify  the  generators  of  conserved  quantities  in 
terms  of  T+  and  T  .  However,  T+  and  T  are  related  to  via 
integral  equations,  and  the  argument  used  in  the  D-bar  approach 
to  show  that  T^L  and  T are  t ime- independent  and  can  be 
expanded  in  powers  of  1/k*  >  can  *lso  be  made  for  the  operator 
version  of  the  problem. 

Quantum  NLS  is  a  1  +  ld  nonrelativistic  theory  for 

complex-scalar  bosons  with  a  four-field  interaction,  or 

equivalently,  rewritten  in  first-quantized  form,  it  describes  a 

system  of  bosons  interacting  via  a  two-body  de 1 ta - f unct i on 

potential.  The  exact  energy  spectrum  can  be  calculated  by  use  of 

Bethe's  Ansatz,  and  this  fact  is  intimately  connected  to  the 
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existence  of  a  corresponding  YB  algebra  '  In  a  sense,  quantum 
NLS  serves  as  a  paradigm  for  other  integrable  quantum  field 
theories  in  1+1  d,  such  as  the  Thirring  model.  We  believe  that 
the  same  will  be  true  for  its  2+ld  generalization, 
Davey-Stewartson,  and  that  the  generalized  YB  relations  will  be 
useful  in  developing  a  generalization  of  Bethes  Ansatz  to  extract 
the  physical  energy  spectrum  of  such  models,  and  also  exact 


S -matr ices . 
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